Electronic structure of LaBr 3 from quasi-particle self-consistent GW calculations 
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Rare-earth based scintillators in general and lanthanum bromide (LaBr3) in particular represent a 
challenging class of materials due to pronounced spin-orbit coupling and subtle interactions between 
d and / states that cannot be reproduced by standard density functional theory (DFT). Here a de- 
tailed investigation of the electronic band structure of LaBr3 using the quasi-particle self-consistent 
GW (QPscGW) method is presented. This parameter-free approach is shown to yield an excellent 
description of the electronic structure of LaBr3. Specifically it is able to reproduce the band gap, 
the correct level ordering and spacing of the 4/ and 5d states, as well as the spin-orbit splitting 
of La-derived states. The QPscGW results are subsequently used to benchmark several computa- 
tionally less demanding techniques including DFT+U, hybrid exchange-correlation functionals, and 
the Go Wo method. Spin-orbit coupling is included self-consistently at each QPscGW iteration and 
maximally localized Wannier functions are used to interpolate quasi-particle energies. The QPscGW 
results provide an excellent starting point for investigating the electronic structure of excited states, 
charge self-trapping, and activator ions in LaBr3 and related materials. 

PACS numbers: 29.40.Mc 71.20.-b 71.20.Eh 71.70.Ej 



I. INTRODUCTION 

Scintillators are materials that exhibit luminescence 
upon excitation by ionizing radiation^ which means that 
a fraction of the absorbed energy is re-emitted as light. 
The emitted photons can be subsequently converted into 
an electric current using for example photomultiplicr 
tubes or photo diodes, which allows one to measure 
the energy spectrum of the incoming radiation. Scin- 
tillation is observed in crystals, plastics, liquids, and 
glasses.— Examples of inorganic crystal scintillators in- 
clude halides, oxides, and chalcogcnidesJ^ The emitted 
light can be produced by exciton recombination (e.g., al- 
kali halides, Csl, MgWCU), core-to- valence, also known 
as cross- luminescence or Auger- free luminescence, (e.g., 
BaF2, CsCl) or most commonly by the relaxation of an 
excited activator atom (e.g., NaI:Tl, Sr±2:Eu, LaF^-.Ce)^- 

Over the course of the last couple years, interest in 
scintillator materials has surged thanks to large scale ap- 
plications in nuclear and radiological surveillance, high- 
energy physics and medical imaging. While the gen- 
eral potential of scintillators has been demonstrated, one 
of the current goals is to develop materials with im- 
proved energy resolution sufficient to detect fissile ma- 
terials with a low probability of errors at ports, borders, 
and airports^ The current state-of-the-art material is Ce- 
doped LaBr3^ which has been extensively characterized 
both experimentally^— and theoretically^— Yet funda- 
mental features of its electronic structure are still incom- 
pletely described and quantified. Such information is, 
however, crucial for understanding the much improved 
performance of LaBr3 compared to other scintillators. 

An ab-initio description of the electronic structure of 
LaBi'3 is challenging due to the presence of La-4/ and 5d 
states as well as pronounced spin-orbit coupling (SOC). 
If one furthermore aims to model Ce activator ions, the 
Ce-derived / and d levels need to be considered as well. 



Previous electronic structure calculations of LaBi-3 
have been based on the Hartree-Fock (HF) method^, 
the LDA+{7 approach^ or hybrid exchange-correlation 
(XC) functionals?2, As is well-known, the HF method 
grossly overestimates the band gap in extended systems 
and the two latter methods both rely on additional fit- 
ting parameters. They are therefore not predictive and 
require experimental information or higher-level calcu- 
lations as reference. In view of this situation the ob- 
jective of the present work is to calculate the electronic 
structure, specifically the quasi-particle (QP) spectrum, 
of lanthanum bromide from an essentially parameter-free 
approach. In doing so one obtains a very good starting 
point for investigating e.g., charge self-trapping, exciton 
spectra or activator level alignment. To accomplish this 
goal the quasi-particle self-consistent GW (QPscGW) 
method^ - — is employed in conjunction with maximally 
localized Wannier functions! 17 ! 18 

The remainder of this paper is organized as follows. 
In the next section the methodology and computa- 
tional details are reviewed. Section MI Al presents as 
the main result of the present work the band structure 
and density of states from QPscGW calculations includ- 
ing SOC. The convergence of our calculations is elab- 
orated in Sect. MI Bl A detailed comparison with both 
hybrid, DFT+J7, and GqWq calculations is presented in 
Sect. MI CI Finally, the main conclusions and an outlook 
of future work are given in Sect. IIVI 



II. METHODOLOGY 

A. The GW approximation 

In the most common implementation of the GW 
approximation^— QP energies are computed to zcroth 



2 



order in perturbation theory according to 
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Here, T, V n - e , and Vh denote the kinetic energy term 
as well as the nucleus-ion and Hartree potentials, respec- 
tively. The renormalization factor (or QP weight) Z n k is 
obtained from the energy derivative of the self-energy.— 
Typically, the unperturbed single particle energies and 
wave functions are obtained from a density functional 
theory (DFT) or HF calculation. This approach is usu- 
ally referred to as the GoWq method. 

Obviously, the application of Eq. ([T]) does not alter the 
underlying single particle orbitals. It is primarily for this 
reason that a starting guess is usually required to yield 
qualitatively the correct level ordering (at least of the oc- 
cupied states). This condition is violated for example in 
lanthanide oxides and can be overcome by generating the 
initial wave function and QP energies using the LDA+C7 
method^ From this description it is apparent that the 
GqWo method yields results that are dependent on the 
initial wave function. 

To overcome these limitations Kotani, Schilfgaardc, 
and Falcev developed the QP self-consistent GW 
method*^— The key idea in the QPscGW approach 
is to optimize an effective self-consistent non-interacting 
Hamiltonian H° = T + V e s such as to reproduce the en- 
ergy dependence of the self-energy as closely as possible 
within the random phase approximation (RPA). This is 
accomplished by minimizing a norm that measures the 
difference between H° and H(uj) with respect to the ef- 
fective one-body potential V e s. A practical scheme is 
obtained by requiring the XC potential to depend on an 
average of the Hermitian part of the self-energy operator, 

V xc = \ { Re P( e i)]« + Re > ( 2 ) 

u 

where Brillouin zone indices have been dropped for 
brevity. It can be argued that within the limits of the 
RPA the solutions of H° can be interpreted as quasi- 
particlesJ£ 

The original form of the QPscGW method was later 
slightly modified by Shishkin, Marsman, and Kresse 2 ^ 
and extended to account for vertex correction in W. For 
practically all materials considered so far the QPscGW 
has been found to give band gaps in very good agree- 
ment with experiment safe for a slight overestimation in 
particular for small gap materials ! 2627 

From our point of view the two major features of the 
QPscGW method in comparison to the GqWq method de- 
scribed by Eq. ([1]) are that (z) the single particle orbitals 
are updated during the course of the self-consistency 
loop and (ii) the final results are independent of the 
initial wave function and thus also independent of any 
adjustable parameters. For the latter statement to be 
fulfilled, it is essential that the basis in which H° is ex- 
panded is sufficiently large to be considered complete as 
discussed in Sect. MI Bl 



B. Maximally localized Wannier functions 

Wannier functions 2 ^ (WFs) provide a complementary 
basis set to the Bloch functions. WFs are defined as 
Fourier transformations of Bloch functions with respect 
to crystal momentum vectors. This formulation is, how- 
ever, arbitrary due to the undetermined phases of the 
Bloch functions. This was utilized by Marzari and Van- 
derbilt who introduced generalized Wannier functions^ 7 - 
for composite bands defined by 
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Here the unitary matrices mix Bloch states having 
the same wave vector k. The quadratic spread of the set 
of WFs can then be expressed in terms of matrix elements 
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and subsequently minimized to obtain so-called maxi- 
mally localized Wannier functions (MLWFs). This for- 
malism was later extended by Souza et al. to the case of 
entangled energy bands.— For the purpose of the present 
work, the underlying first-principles code was modified to 
enable the calculation of the matrix elements appearing 
in Eq. which were then used as input for wan- 

nier9022 to obtain MLWFs. The site and angular mo- 
mentum projected QP energies were finally interpolated 
using the MLFW basis to generate accurate band struc- 
tures and DOS. 



C. Spin-orbit coupling and symmetry operations 

The spin-orbit correction to the Hamiltonian can be 
obtained from the Dirac equation by three successive 
Fouldy-Wouthuysen transformations^ In atomic units, 
the resulting Hamiltonian for an electron in a central po- 
tential V can be written as 



2 r dr 



(5) 



Here a is the fine structure constant and t and s 
are the orbital and spin angular momentum operators, 
respectively. Aryasetiawan and Biermann generalized 
Hedin's equation s) 19 ' 20 to the case of spin-dependent 
interactions ! 32 ' 33 The algebraic structure of the modi- 
fied set of equations is very similar to Hedin's original 
equations, although it is shown that in general the self- 
energy becomes spin-dependent and the polarization now 
describes the response of the charge density to magnetic 
fluctuations and vice versa. In the absence of explicit 
two-particle spin-interaction the self-energy will only de- 
pend on variations in the electric field.— Therefore, if 
vertex corrections, spin-spin, and spin-other-orbit inter- 
actions are neglected, this reduces for a non-magnetic 
system to adding H so to the unperturbed single-particle 
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Hamiltonian and embodies an explicit spin dependence 
of the Green's function. Sakuma et alM' applied this for- 
malism to study the effect of SOC in Hg chalcogenides in 
the Go Wo approximation and found an enhancement of 
the spin-orbit splitting of 0.1 eV. In the context of QP- 
scGW calculations Chantis et al. have included SOC as a 
perturbation on top of already converged orbitals for III- 
V and II- VI zincblende semiconductors^ In the present 
case, to investigate spin-orbit splitting SOC is included 
in the self- consistency cycles of the QPscGW method. 

In the projector augmented wave metho d 35 ' 36 (PAW) 
the SOC Hamiltonian consists of three terms 

H so = H so + ^ \pi) ((4>i\H so \(f)j) - ((f>i\H so \(t>j)) 

ij 

(6) 

where the projectors p and orbitals <f> and <p have their 
usual meaning^ In DFT calculations it is customary to 
only include the SOC in a sphere around each atom, us- 
ing the spherically averaged self-consistent Kohn-Sham 
potential. Assuming a complete PAW basis set, the first 
and third terms in Eq. © cancel. In the current imple- 
mentation only the second term is kepti 3 -^ 

For the QPscGW calculations reported in this work 
the local potential in Eq. ([5]) is computed at each step in 
the self-consistent cycle as the spherically averaged sum 
of the external, Hartree and PBEiS XC potential^ As 
a result, these calculations include SOC only on a DFT 
level. 

To render hybrid XC and GW calculations including 
SOC computationally feasible it is imperative to take ad- 
vantage of the space group symmetries and solve for QPs 
within the irreducible Brillouin zone only. To this end we 
augmented the action of a symmetry operator in Carte- 
sian space on wave functions, projected wave functions, 
and the action of the self-energy operator and its deriva- 
tive with respect to energy on a Bloch function to include 
rotations in spin space. Now, a general space group oper- 
ation contains a point group operator S and a translation 
w. The action of this operator on a spinless Bloch func- 
tion isM 



are the eigenvector matrix and diagonal eigenvalue ma- 
trix of A, respectively. 

We also note in passing that operations of the type 
given by Eq. ([7} coupled with spin rotations are also use- 
ful for the evaluation of the matrix Af( k ' b ) defined in 
Eq. Q. 

D. Computational details 

Calculations have been performed using the projector 
augmented wave metho d 35 ' 36 as implemented in the Vi- 
enna ab-initio simulation packaged - — First wave func- 
tions and QP energies were obtained within a general- 
ized Kohn-Sham scheme using several different approx- 
imations to represent exchange and correlation effects. 
Subsequently the single particle states of these calcula- 
tions served as starting points for Go Wo and QPscGW 
calculations. For the XC potential in the initial calcula- 
tions, we considered a generalized gradient functional^ 
the BFT+U method^ with parameters from Ref. ka 
several range-separated hybrid XC functionals with vari- 
able mixing parameter a and a fixed screening length of 
fx = 0.2 A |2£ as well as exact-exchange (EXX), where 
the latter is equivalent to carrying out a restricted HF 
calculation. 

All calculations were performed at the experimental 
lattice parameters^ using a T-centered 6x6x3 grid, a 
general plane-wave cutoff energy of 219 eV, and a cut- 
off energy of 146 eV for the response function in the 
GW loop. The Green's function and the screened in- 
teraction in the GW loop were evaluated using 2088 and 
1056 bands in calculations with and without SOC, re- 
spectively. The effective QPscGW Hamiltonian was ex- 
panded in a basis containing up to 192 bands (384 when 
including SOC) equivalent to states up to 34 eV above 
the conduction band minimum (CBM). In general, our 
parameters ensure convergence of the valence band and 
lower conduction band QP energies of about 0.05 eV. The 
convergence of our calculations is demonstrated and dis- 
cussed in Sect. IIII Bl 



{S'liu} exp [iSk ■ r] Uf.(r) = 

exp[iSk ■ (r — w)]usk(r — w), (7) 

where Uu, is the cell-periodic wave function. The rotated 
wave function thus corresponds to a Bloch function of 
wave vector Sk. If the point group operator S contains 
a rotation by an angle a around some unit axis n the 
action of the space group operator on a Bloch spinor 
would be augmented by an additional rotation of the two 
spin components by the spin-space rotation matrix 



R s (n, a) = exp 



ictcr ■ n 



exp(A) 
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where <7j is a Pauli matrix. In practice this rotation ma- 
trix is computed by R s = V' cxp(D)V where V and D 



III. RESULTS 

A. Band structure and density of states 

Figure [1] represents the main result of the present pa- 
per. It shows the band structure and density of states 
(DOS) of lanthanum bromide from QPscGW calculations 
taking SOC into account. As in virtually all halides that 
we are aware of, the valence band is predominantly com- 
posed of halogen p-states. Similar to other compounds 
with comparably low symmetry and a valence band max- 
imum (VBM) that is prim arily derived from p-states (for 
example In 2 03, Refs. l49l andlHol. and Sri2, Ref. [sT|) , the 
VBM in LaBr3 is very flat indicating a very low hole 
mobility. 
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FIG. 1. (a) Band structure and (b) density of states of LaBr3 
as obtained from QPscGW calculations including spin-orbit 
coupling. Experimental X-ray photoelectron spectrum from 
Ref . EE 
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FIG. 2. Comparison of (a) valence and (b) conduction band 
density of states from QPscGW calculations with and with- 
out spin-orbit coupling (SOC). Valence and conduction band 
states are shown with respect to valence band maximum and 
conduction band minimum, respectively obtained from calcu- 
lations without SOC to highlight the SOC-induced shift. 



By comparison the conduction band structure is more 
complex. Its bottom mostly consists of La-5<i states that 
extend up to about 3.0 eV above the CBM where a mini- 
mum in the DOS is observed. At slightly higher energies 
a pronounced peak due to La-4/ states is clearly visible, 
superimposed on a rather broad band that has predomi- 
nantly Br-derived d-character. 

The SOC lifts the top of the valence band by 0.19 eV 
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FIG. 3. (a) Absorption spectra measured by Sate— - and (b) 
state resolved contributions to joint density of states (JDOS) 
from QPscGW calculations due to transitions from valence 
band (VB) states to different conduction states (SOC contri- 
butions are included). 



while leaving both the bottom of the Br-4p band and 
the Br-4s band unchanged. The most significant effect 
is observed for the La-4p derived states which split by 
2.25 eV in excellent agreement with experiment [com- 
pare Fig. [Ha)]. In general over the energy considered 
the agreement of the QPscGW calculation with X-ray 
photoelectron spectroscopy (XPS) data^ is very good if 
SOC is included. 

The bottom of the conduction band, which is com- 
posed of La-5d states, is virtually unaffected by SOC. 
There is also no change in the position of the La-4/ peak 
situated about 3.2 eV above the CBM, although one no- 
tices a slight redistribution of weights in this band. The 
QPscGW calculations yield a band gap of 6.19 eV, which 
is reduced to 5.99 eV when SOC is taken into account. 
This value is in very good agreement with experiments, 
which indicated a value of 5.9 eV£ 

To the best of our knowledge, there are currently no 
spectroscopic data corresponding to the single-particle 
excitation spectrum that would enable a direct compar- 
ison of the calculated conduction band structure with 
experiment. Sato^ measured the absorption spectrum, 
which depends on both valence and conduction band 
states. Since excitonic effects are known to be impor- 
tant in halides, as shown specifically for LaBr3 in Ref. [f| 
an accurate description of absorption would require the 
inclusion of electron-hole interactions (see e.g., Refs. l52l 
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and[53|). This could be accomplished for example by solv- 
ing the Bethe-Salpeter equation^ - — which is, however, 
the subject of future work. In the present paper, we re- 
sort to a simplified comparison based on the joint density 
of states (JDOS). On a single-particle level, the latter is 
related to the dielectric function and thus the absorp- 
tion coefficient as discussed for example, in Ref. l58l . One 
can therefore expect a correlation between characteristic 
features in the absorption spectrum with features in the 
JDOS. 

As observed in Fig. [31 there indeed exist similar fea- 
tures in absorption spectra measured by Sato^ and the 
JDOS obtained from QPscGW calculations, as indicated 
by letters A and B. 

When separating the contributions from different con- 
duction band states to the JDOS, FigureEJb) shows that 
features A and B arise from transitions from valence band 
(VB) to La-d and / states, respectively. The fact that the 
computed features occur at slightly higher energies com- 
pared to experiment is consistent with the observation 
that electron-hole interactions usually lead to a red shift 
of the calculated spectrum. In Sect. IIII CI we compare 
the JDOS obtained from different computational meth- 
ods, which will demonstrate that reproducing the two 
features discussed above is not trivial. The good agree- 
ment between the QPscGW JDOS and the XPS data 
therefore provides evidence for the reliability and accu- 
racy of the QPscGW results. 



B. Convergence of QPscGW calculations 

In Sect. |n] it was argued that the QPscGW method 
yields results that arc independent of the starting wave 
function. This feature is illustrated in Fig. Uta), which 
compares the convergence of several characteristic QP en- 
ergy differences starting from both DFT and EXX wave 
functions. The two different starting points correspond 
to a significant underestimation (DFT) and overestima- 
tion (EXX) of the band gap, respectively. Yet after con- 
verging the QPscGW cycles, the converged QP energy 
differences agree to within 0.05 eV for valence and lower 
conduction band levels, and to 0.18 eV for the La-4/ 
states. The remaining differences can in principle be re- 
duced further by increasing the basis set in which the ef- 
fective Hamiltonian of the QPscGW method is expanded. 

The previous statement is confounded by Fig. HJb), 
which shows the convergence of several QP energy differ- 
ences with respect to the number of bands included in the 
QPscGW basis set £2 Due to the presence of both La and 
Br d-states as well La-/ states, the density of states in 
the conduction band is larger than in conventional semi- 
conductors such as Si or GaAs. As a result, one needs five 
times as many unoccupied than occupied states in order 
to reach the convergence exhibited in Fig.[4ta), a number 
that is considerably larger than for the aforementioned 
semiconductors. In general we have found it very useful 
not only to check convergence with respect to the num- 
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FIG. 4. (a) Convergence of QPscGW calculations start- 
ing from DFT and EXX wave functions, respectively, as evi- 
denced by several quasi-particle energy differences, (b) Con- 
vergence with respect to QPscGW basis set size starting from 
DFT wave function. The results shown in (a) were obtained 
using 192 bands to expand the effective QPscGW Hamilto- 
nian, corresponding to states up to 34 eV above the CBM. 
The calculations shown do not include SOC. 



ber of bands as shown in Fig. |U(b) but also to confirm 
convergence of our calculations by comparing the results 
obtained from initially very different wave functions such 
as DFT and EXX. 



C. Hybrid DFT and G Wo calculations 

For studying various problems of interest such as 
charge self-trapping, defect formation, or alignment of 
activator levels one needs to consider ionic relaxations 
and representative supercells. At present QPscGW cal- 
culations are, however, computationally still extremely 
expensive because a large number of unoccupied bands 
has to be included. More severely, the method does not 
allow to obtain total energies and forces. It is therefore 
important to determine to which extent computationally 
lesser demanding calculation schemes can reproduce the 
reference electronic structure provided by QPscGW cal- 
culations. To this end, a number of conventional and 
hybrid exchange-correlation functionals, as well as the 
corrections from Go Wo to some of these functionals have 
been considered. 
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FIG. 5. (a) Band gap, (b) position of La-4/ states, and (c) 
valence band width from conventional exchange-correlation 
functional (PBE), several hybrid exchange-correlation func- 
tional (hybrid XC), and exact-exchange (EXX) calculations. 
Vertical arrows indicate the improvement obtained by carry- 
ing out Go Wo calculations. Gray bars indicate QPscGW val- 
ues that serve as reference data in this comparison. All data 
obtained without spin-orbit coupling. In (a) two values are 
shown for PBE corresponding to the gaps between VBM and 
La-4/ and 5d states, respectively (compare Fig. [6]). For all 
other cases the band gap between VBM and La-5d is shown. 



Figure [5] shows the variation of several QP energy 
differences with the mixing parameter a of a range- 
separated hybrid exchange-correlation functional^ with 
a screening length of /Lt = 0.2 A -1 . For a = this 
functional is identical to the XC functional by Perdew et 
al. (PBE)i^ For comparison we also included data from 
EXX calculations, which do not include correlation. 

The data shows that at the PBE level, the ordering of 
the La-4/ and 5d states is incorrect with 4/ located be- 
low 5d levels. Upon inclusion of exact-exchange and/or 
by adding a Go Wo calculation the level ordering is cor- 
rected but for most functionals the QP energies are still 
quantitatively very different from the QPscGW values. 
In general as a is increased the QP energy differences con- 
sidered here increase. Performing Go Wo calculations in 
general improves the agreement. The same trend was ob- 
served by Aulbur et al. for medium and wide-gap materi- 
als using Go Wo on top of an EXX-based XC functional^ 
While the Go Wo method reduces the dependence of the 
results on a, it remains pronounced. 
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FIG. 6. Density of states as obtained from various different 
calculation schemes (without SOC). 



For a values between approximately 0.5 and 0.65 all 
QP energy differences considered in Fig. [5] intersect the 
respective reference values. This suggests that a hybrid 
functional with a in this range could possibly reproduce 
the DOS obtained from QPscGW calculations. We there- 
fore considered a hybrid functional with a = 0.62, for 
which the total DOS is shown in Fig. [5] together with re- 
sults obtained from PBE, EXX, DFT+(7, and QPscGW. 
For simplicity SOC effects were neglected in this com- 
parison. In fact apart from the deep semi-core states the 
DOS for the modified hybrid functional compares very 
favorably with the QPscGW reference. 

The figure also contains the DOS obtained from 
DFT+U calculations using U and J parameters deter- 
mined previously for other La-compounds^ both as- 
calculated (DFT+U) and with a rigid shift of the con- 
duction bands to the QPscGW band gap ( "scissors" cor- 
rection, DFT+U + x). Both the upper valence and lower 
conduction band structures from DFT+U + \ are in ex- 
cellent agreement with QPscGW data, and even though 
the agreement worsens for deeper lying valence states, 
among the methods considered the DFT+U +\ approach 
still yields the best agreement with the reference DOS. 
Still for practical calculations one would have to resort 
to the DFT+U method without the scissors correction, 
which exhibits a pronounced band gap underestimation 
(3.60eV vs 6.19eV without SOC). 

In this context it is furthermore interesting to revisit 
the JDOS, shown in Fig. [7J which was compared earlier 
with absorption data (see Fig. [3]). Also in this arena both 
the DFT+U method and the modified hybrid XC func- 
tional (not shown in Fig. [7J yield good agreement with 
QPscGW calculations exhibiting the features discussed 
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FIG. 7. Joint density of states from different calculation 
schemes (without SOC). It is instructive to compare these 
curves with the data in Fig. [3] 

in detail in Sect. IIII Al In stark contrast the JDOS from 
PBE and EXX calculations are qualitatively very differ- 
ent from the QPscGW data. The position of the features 
in the absorption spectra shown in Fig. [3] are not repro- 
duced with these functionals. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper it has been demonstrated that QPscGW 
calculations of LaBr3 yield single particle spectra (DOS) 
that closely match experimental XPS data (Fig. [T]). At 
the same time the key features in the JDOS correlate with 
those in the experimental absorption spectrum (Fig. [3J . 
To obtain this result it was essential to include spin- 
orbit coupling in the calculations. The very good agree- 
ment between our QPscGW calculations and experimen- 
tal data provides confidence in the reliability of the QP- 
scGW method for treating systems with /-electrons and 
weak to moderate correlation. Our work complements 
earlier studies of other /-electron systems that used both 
the QPscGW methodic and the G W approach^ but 
focused on systems, in which /-electrons are located in 
the valence band. 

The presence of both / and d states in the conduction 
band of LaBr3 requires — by comparison with more con- 
ventional semiconductors and insulators — an unusually 
large basis set for the expansion of the effective Hamilto- 
nian to obtain converged results [Fig. H£b)]. It was also 
explicitly demonstrated that the QPscGW results are in- 
dependent from the initial wave function [Fig.[4ta)]. This 
property is a big advantage of the QPscGW approach 
compared to other GW methods, as exemplified in Fig.[5j 
which shows the Go Wo results to exhibit a rather pro- 
nounced dependence on the initial wave function. This 
behavior is important to keep in mind when interpreting 



such calculations. 

Another aspect that deserves mentioning is related to 
the computation of matrix elements and their interpola- 
tion. The QPscGW method does provide updated wave 
functions that in turn can be used to compute matrix 
elements. In contrast the wave functions are unchanged 
when doing Go Wo calculations. This difference becomes 
particularly apparent when interpolation methods are 
used. These methods, based for example in the present 
work on maximally localized Wannier functions, can be 
used to represent the 1 aforementioned matrix elements on 
a finer mesh in reciprocal space, dramatically improv- 
ing convergence of the calculations. If, as in the case of 
LaBr 3 , the level ordering changes as the result of a Go Wo 
calculation the connectivity of the states changes as well, 
which causes problems for interpolation methods. In such 
a case the interpolation of the QP energies and the pro- 
jections on the basis of a PBE + Go Wo calculation are 
physically not meaningful because the derivative matrices 
7y/( fc ' b ) n0 longer refer to the same states. The QPscGW 
method, however, does not suffer from this shortcoming 
and matrix elements can be readily interpolated once the 
MLWFs have been determined. 

It is instructive to compare the electronic structure of 
LaBi - 3 both with other lanthanum halides and the free 
La atom. As discussed at length in this paper in LaBi-3 
the La-4/ levels are located above the La-5cZ levels. This 
is similar to the ordering of excited states in a free La 2+ 
ion£2- In contrast for the 3+ charge state of the La ion 
the excited 4/ states are indeed observed to lie below 
the 5d levels. This situation is compatible with a partial 
charge transfer between La and Br. 

Finally, it was shown that a hybrid XC functional can 
be constructed that yields reasonable agreement with 
DOS (Fig. El) and JDOS (Fig. [7jj obtained from QP- 
scGW calculations. This functional should be suitable 
for studies that require ionic relaxation and/or the use 
of supercells. The DFT+t7 method (using parameters 
from the literature^) also reproduces the DOS structure 
of valence and conduction band states independently but 
still underestimates the band gap significantly. 
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